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The quasi-geostrophic two-layer model is of superior interest in dynamic meteorology since 
it is one of the easiest ways to study baroclinic processes in geophysical fluid dynamics. The 
complete set of point symmetries of the two-layer equations is determined. An optimal set 
of one- and two-dimensional incquivalent subalgcbras of the maximal Lie invariance algebra 
is constructed. On the basis of these subalgcbras we exhaustively carry out group-invariant 
reduction and compute various classes of exact solutions. Where possible, reference to the 
physical meaning of the exact solutions is given. In particular, the well-known baroclinic 
Rossby wave solutions in the two-layer model are rediscovered. 

1 Introduction 

There is a long history in dynamic meteorology to use simplified models of the atmosphere rather 
than the complete set of hydro-thermodynamical equations to study only selected phenomenon 
instead of accounting for the whole variety of weather and climate at once. The greatest simpli- 
fication which is still capable to qualitatively (and, under some conditions, also quantitatively) 
describe the behavior of large-scale geophysical fluid dynamics is the barotropic vorticity equa- 
tion. While this equation indeed allows one to explain the propagation of the mid-latitude 
Rossby waves, it cannot be used to elucidate the occurrence of developing weather regimes. 
The reason for this substantial lack is that the barotropic vorticity equation is a single equa- 
tion valid only in one atmospheric (pressure) layer. However, development in the atmosphere is 
usually associated with the vertical structure of, e.g., the entropy field and hence a single-layer 
consideration is at once limited. 

Of course, the atmosphere is continuously stratified and hence it is in fact three-dimensional. 
However, the main process of baroclinic instability, which is the dominant mechanism responsible 
for the formation of mid-latitude weather systems can already be qualitatively understood by 
considering only two coupled atmospheric layers. Moreover, the studies of layer models create a 
basis for the more complicated investigation of the three-dimensional governing equations. For 
this reason, many results for the two-layer model are available in dynamic meteorology [13, 23], 
making this model particularly interesting for a systematic mathematical investigation. 

The aim of this paper is to carry out Lie symmetry analysis of the quasi-geostrophic two-layer 
model. There already exist a number of papers dealing with symmetries and exact solutions 
of the simpler barotropic vorticity equations [4, 5, 15, 16, 17] as well as the more complicated 
Euler or Navier-Stokes equations (see, e.g., [2, 9, 19, 20, 26] and references therein). At the 
same time, the intermediate two-layer model has not been considered in this light so far. 

The organization of the paper is the following: The model equations are presented in Sec- 
tion 2, together with some of their known properties. Section 3 contains the results on Lie 
symmetries of these equations. The theorem on the complete point symmetry group of the two- 
layer model is proved in Section 4. A list of one- and two-dimensional inequivalent subalgebras 
of the corresponding maximal Lie invariance algebra is constructed in Sections 5. In Sections 6 
and 7, systems obtained under reduction upon using the constructed one- and two-dimensional 
subalgebras are derived and investigated. The Rossby wave solutions in the quasi-geostrophic 



two-layer model are recovered as Lie invariant solutions. Related boundary- value problems are 
discussed in Section 8. Section 9 summarizes the most important results of this paper. Finally, 
in the Appendix A we give a new symmetry interpretation of a method for finding exact solutions 
of linear systems of PDEs. 

2 The two-layer model 

The first impulse to the investigation of baroclinic instability in the two-layer model was given 
in [24]. Since two layers are considered, the model is capable of studying the interaction of 
the barotropic mode and the first baroclinic mode, which is sufficient for describing the basic 
mechanism of baroclinic instability. The model consists of two copies of the barotropic potential 
vorticity, evaluated at two different atmospheric levels [23]: 

Ql + {^\Q l } = o, 
Q 2 t+{*P 2 ,Q 2 } = o, 

where ip 1 and tp 2 are the stream functions in the upper and lower layer, 

Q^vV + zfy-i^ 1 - A 

Q 2 = \J 2 ^ 2 + /3y + F^j 1 - ip 2 ), 

are the corresponding potential vorticities, with the constants f3 and F being the Rossby pa- 
rameter and internal rotational Froude number, respectively, and 

W,Q'} = i>:Ql-i>lQl 

is the usual Poisson bracket of the functions ip' and Q* with respect to the space variables x 
and y. For simplicity, we set the Froude numbers of the two layers to be equal, i.e. F\ = F2 = F, 
thereby assuming both layers to be of equal depth. Moreover, flat topography is assumed. For 
the configuration to be stably stratified, the lower layer must be denser than the upper layer. 

Here and in what follows subscripts of functions denote differentiation with respect to the 
corresponding variables and prime denotes differentiation with respect to t. 

Due to equivalence transformations of scaling and alternating signs in the class of equations 
of form (1), it would be possible to set F 6 { — 1, 0, 1} and j3 £ {0, 1}. Since F and j3 are positive 
meteorological quantities, it would imply that F = ft = 1 but we will not use this scaling below. 

3 Lie symmetries 

In the case F/3 / 0, which is the single subject of the present paper, system (1) admits the 
maximal Lie invariance algebra q generated by the following basis elements: 

d t , d y , X{f) = f(t)d x -f'(t)y (8^,1+8^2), T = d^i-d^, Z(g) = g(t)(d^i +d^t), (2) 

where / and g run through the set of real- valued functions of t. The computation of the algebra q 
was checked using the symbolic calculation packages DESOLV [8] and muLie [11]. For g to really 
be a Lie algebra, additional restrictions on the smoothness of the parameter functions / and g 
should be imposed [9]. 

Physically, these generators are exponentiated to give time and north-south translation sym- 
metry, generalized Galilean transformations with respect to x, as well as shifts and gauging of 
the stream functions. 

The structure of the Lie symmetry algebra suggests the introduction of the new dependent 
variables ip + = ip 1 -\- ip 2 and ip~ = ip 1 — ip 2 . This change transforms the set of generators to 

dt, d y , X{f) = f(t)8 x -2f'(t)yd^ + , T = 1d^-, Z(g) = 2g(t)8^ + . 
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From the meteorological point of view, the new variables have a sound physical meaning. Since 
the two-dimensional wind fields at both the levels can be represented as the derivatives of the 
respective stream functions, it follows that the component ip + gives rise to the mean of these 
fields. This part is usually referred to as the barotropic part of the flow. In turn, derivatives 
of tl>~ give rise to the difference in the wind field between the two-layers. In meteorology, this 
difference is called the thermal wind, which is a measure of the baroclinity of the fluid. Therefore, 
these variables are commonly used in the investigation of baroclinic instability, e.g., in the study 
of the linearized two-layer model [13]. However, from the group-theoretical point of view, their 
usage is already suggested by the special form of Lie symmetry operators (2). 
Using the variables ip + and tp~ , the model (1) is represented as 

Qt + l({^ + ,Q + } + {^~,Q~}) =o, 

7 (3a) 

Qt+ 2 {{^ + ,Q~} + {^~,Q + }) =o> 

where 

Q+ = VV + 2/3y, 

9 (3b) 

are the barotropic and baroclinic potential vorticities, respectively. In the present paper, we 
will use both the forms of the two-layer model, i.e. employing both the "layered variables" if) 
and ijj 2 and the "barotropic/baroclinic variables" tp + and ip~ . The precise usage depends on 
whether we study linear or nonlinear submodels of the two-layer model. 



4 Complete point symmetry group 

The complete point symmetry group of a system of differential equations, which includes both 
continuous and discrete symmetries, is conventionally calculated using the direct method. The 
outlines of this method are quite simple. Supposing the most general form of a point transforma- 
tion in the associated space of independent and dependent variables, one expresses all involved 
derivatives of the new (transformed) dependent variables via the old variables, substitutes the 
obtained expressions into the system written in terms of the transformed variables, then ex- 
cludes the derivatives which are assumed constrained due to the system (principal derivatives) 
and splits with respect to the unconstrained (parametric) derivatives. As a result, one obtains a 
system of determining equations for point symmetry transformations, which is nonlinear in con- 
trast to a similar system arising under application of the infinitesimal Lie method and, therefore, 
is much more complicated for solving. This is why different special techniques (the implicit rep- 
resentation for unknown functions, the combined splitting with respect to old and new variables, 
inverse expression of old derivative via new ones, a mapping of the initial system by a point 
transformation to another system, etc.) are applied to the derivation of determining equations 
and their a priori simplification [6, 27, 30, 31]. 

Here we aim to use an approach similar as described in [14], which is based on the knowledge 
of the Lie symmetries of a given differential equation. This method rests on the fact that any 
point symmetry generates an automorphism of the maximal Lie invariance algebra. By factoring 
out the continuous symmetries from the whole point symmetry group, the discrete symmetries 
of the differential equation can be determined. 

The computation of the complete point symmetry group can be considerably simplified by 
noting that any automorphism of a Lie algebra g leaves invariant, by definition, all megaideals 
of g. Recall that a megaideal of g is a vector subspace of g which is invariant under any transfor- 
mation from the group of automorphisms of g [28]. Therefore, by determining megaideals of the 
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maximal Lie invariance algebra of the given differential equation and imposing the invariance 
condition of these megaideals under the push-forwards of vector fields associated with the point 
symmetries allows to restrict the general form of possible point symmetries already before trans- 
forming the differential equation itself. After taking into account these initial restrictions, it is 
usually much simpler to proceed with the splitting of the variables in the transformed differential 
equation as described in the first paragraph. 

In this section, the approach just outlined is demonstrated for the two-layer equation in 
barotropic/baroclinic variables. 

Theorem 1. The point symmetry group G of system (3) consists of transformations of the form 

t = eit + T , x = eix + f(t), y = e 2 y + Y , 

#~ = + , i> + = e 2 ?p + - lE X e 2 fty + g(t), 

where £j = ±1, i = 1, 2, 3; Tq, Yq, £ 1R and f and g are arbitrary smooth functions of t. 

Proof. Recall that the maximal Lie invariance algebra of system (3) is the infinite dimensional 
algebra g = (dt, d y , X(f),J-, 2(g)), where / and g run through the set of smooth functions of t. 
It is a solvable algebra since q' = (X (/), 2(g)) and hence q" = {0}. In other words, the radical r 
of g coincides with the entire g. The algebra g has the nontrivial center 3 = (X(l), J 7 , 2(1)). 

The nil-radical of g is the ideal n = (d y , X(f), J 7 , 2(g)) . Indeed, this ideal is a nilpotent 
subalgebra of g since n^ 2 ) = n' = [n, n] = (2(g)) and n^ 3 ) = [n, n'] = 0. A unique ideal of g 
properly containing n is the entire algebra g itself, which is not nilpotent. This means that n is 
the maximal nilpotent ideal. 

In the calculations of G, we use the following megaideals of g: the entire algebra g, the 
derived algebra g', the nil-radical n, its derivative n', the center 3 and their proper intersections, 
3 n g' = (X(l),2(l)) and 3 n n' = (2(1)). 

For the general point transformation 

T: (t,x,y,4>-,4> + ) = (T,X,Y,^,t> + ), 

where T, X, Y, ^ + are functions of t, x, y, ip + with the Jacobian not equal to zero, 
to be a point symmetry of system (3), the associated push- forward 71 of vector fields must 
be an automorphism of g. In particular, %Q = g, %q' = g', %n = n, %n' = n', %$ = 3, 
%{l n g') = 3 n g' and 7; (3 n n') = 3 n n'. 

Investigating the restrictions on T imposed by the invariance of 3 n n' under T*, we have 

%2(1) = 2(V^~ + X^d, + Y^dy + *- + ^_ + *J + V) = ^( a )» 

where a = const. This equation implies that = X^+ = Y^+ = = 0, ^ + = a = const, 
and Then, from the transformation of the elements of n' we conclude 

%2(g) = 2agd^ + = Z(g% 

where g g is a smooth function of t related to g. Comparing coefficients, we find that ag(t) = 
g 9 (T). As g is an arbitrary smooth function of t, we can set g = t to obtain 2at = g l (T). 
Because of a 7^ and hence g~ 7^ 0, this implies that T = T(t) and Tt 7^ 0. 
In a similar manner, the condition 

% X(l) = X(h)+ Z(b 2 ), 

with b\,b 2 = const follows from the invariance of 3 n g' with respect to T*. It is spit into the 
equations Y x = \E'~ = 0, X x = b\ = const and = 2b 2 = const, where b\ 7^ 0. The invariance 
of g' implies the condition 

%X(f)=X(p) + 2(gf), 
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where y and gf are smooth functions of t related to /. Comparing coefficients in the last 
condition immediately gives the additional equations b±f = ff(T) and b 2 f — aftU = —f({T)Y + 

gf(T). Taking into account the equality f~ = bift/Tt, we obtain that Y = Y 1 (t)y + Y°(t), 
where Y 1 = aT t /b\ 7^ and the precise expression of Y° is not essential for this time. 
The push-forward of the remaining basis operator T of the center 3 by T implies 

%T = X{ Cl ) + c 2 T + Z{c 3 ), 

where again ci, 02,03 = const. From this condition, we can conclude that X^- = c\ = const, 
\E^_ = c-2 = const and \E'^_ = 2C3 = const. 

It remains to investigate the restrictions imposed by the push-forwards of the basis operators 
dt and d y . The operator dt does not lie in the above proper megaideals. Hence, its push-forward 
can be represented only as a general element of g: 

%d t = d x di + d 2 8y + d 3 P + X{f) + Z(g) 

for real constants d\, d 2 , d 3 and smooth functions / and g of t. It is then straightforward to find 
Tt = d\ = const, Yt = d 2 = const and thus Y 1 = const and Y° = d 2 t + d^, where d± = const. 
Moreover, ^ = 2d 3 = const, X t = f(T) is a function of t and = -2f~ t {T)y + g(T). 

Since the operator d y belongs to the megaideal n, its push-forward by T should take the form 

%d y = etdy + e 2 F + X(f y ) + Z{g y ) 

with ei,e2 = const and f y and g y again being smooth functions of t. Comparing coefficients 
implies Y 1 = e 1} ¥~ = 2e 2 , = -2f?{T)y + 2g y (T) and X y = f y (T), which is a function of t. 

As X ty = 0, Xy = const holds. Therefore, f~ = and depends only on t. 

Collecting all constraints obtained so far and re-denoting the involved values, we obtain the 
representation of the point transformations inducing automorphisms of g: 

T = T 1 t + T , X = X lX + X 2 y + X^~ + /(t), Y = Y iy + Y 2 t + Y , 

q>- = + + ^3 t + % , ( 4 ) 

= + + *3 f ^ _ + <p(t)y + g(t), 

where T , T 1; X x , X 2 , X 3 , Y , Y u Y 2 , f 3 ", are constants, T^Y^^ + 0, 

X\Y\ = *ifj T\, f and g are smooth functions of t and ^ = —2Y\f t t- 

We now have to take into account that the transformation T is a point symmetry of sys- 
tem (3). Therefore, we should find explicit expressions for the derivatives of ip~ and ip + with 
respect to the new variables t, x, y. In view of the representation (4), the transformation rules 
for the partial derivative operators read 

1 / D t X Y 2 ( B V X 
d T = — D 4 —B x - — ( D„ v — D, 

^ = — !— a 5 = — (d„ - 
x I),,Y ' ' V, V D, V 

where D^, D x and Dy denote the operators of total differentiation with respect to t, x and y, 
respectively. The derivative 4>txy can ar i se on ly i n the expression for 

7+ ^ 1 + 

^ra = ~TxY^ T\X^ txy + " ' ' 

Since in the first equation of (3) there is no term with ip^. y , we consequently have ^> y X = 0, 
which implies X 2 = X3 = 0. Using this result, the transformation rules for the partial derivative 
operators are essentially simplified: 

&i = 7 k( dt ~ x\ Qx ~ y\ dy ) 1 9& = A^' &i = Y\ Qy ' 
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Upon the substitution ip~£ x = V 2 ip + — ip£ y , we obtain 

Since there is no extra term ipf in the first equation of (3), we have X\ = Y 2 =: a ^ 0. 

To plug the transformed variables into system (3), it is convenient to write down the expres- 
sions for the transformed potential vorticities Q~ and Q + : 

Q + = —(Q + - 2(3y) + ^ (Q- + 2F4>-) + 2p{Y x x + Y 2 t + Y ), 
a a 

QT = -^(Q~ + 2F ^~) - 2F{^ir + ^y + ^3 t + ). 

The further restrictions on the transformation T can be found by splitting the resulting equations 
with respect to the variables t, x, y, t/) x , ipy, Qt an d ■ 

We start with the restrictions imposed from the transformation of the second equation of 
system (3). The term with ip~Q~ arises after the expansion of V'J'Q^ + ^JQ^ ■ As VD^ / 0, 
the corresponding coefficient, which is equal to \Eg /(gX\Y\), vanishes if and only if = 0. 
A term with ip^ is contained only in the expression for . The corresponding coefficient 
2FT 1 _1 ^^((j _1 — 1) also must be equal to zero, i.e., a = 1 and hence X\ = z\ = ±1 and 
Y± = £2 = ±1. Splitting in a similar way with respect to ip^ an d ipy , we respectively find that 
^ 2 =0 and Yi = i &f. As we already know that Y 1 = ^fTi/Xi, this implies that T\ = X\ = e\. 

From these restrictions, we can conclude the following form of the transformations: 

T = e l t + T , X = s lX + f(t), Y = e 2 y + Y 2 t + Y , 

= + q>~t + 9q , = e 2 ip + + $>t x + <p(t)v + 9 (*), 

Q~ = *J"Q- - 2F(^3 t + *q ), Q + = e 2 Q + + 2/3(y 2 t + Y ). 

Substituting these transformations into the second equation of system (3), we find 

^- (QT - f -rQT x - ^Qy) - + J—( £ rt* + *t)*IQy - Mt + <p)*IQ*) 

+ f^* - ^-q+) = |i- (V + - vtf*rQ«) + ^* - ^r<#)) • 

Splitting of this equation immediately gives ^3 = 0, (p = —2e^ 1 e 2 ft and Y 2 = \^ 2 . 

In a similar fashion, plugging these transformations into the first equation of system (3) we 
obtain 

J- (Qt ~ T Q * ~ ^ Q v) + + oTT^* + *^> 2 Q+ - + <P)S2Qt) 

£1 \ £1 £2 / £1 2eie 2 

The symmetry condition implies = £3 = ±1, ^2 — and consequently = 0. This 
completes the proof of the theorem. □ 

Remark 1. The continuous transformations generated by elements of the center 3 and only such 
transformations from the point symmetry group of system (3) induce the identical automorphism 
of the algebra g. 
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Remark 2. Comparing the results presented in Theorem 1 and Section 3 leads to the conclusion 
that besides Lie point symmetries system (3) admits discrete point symmetries. The group of 
discrete symmetries is generated by the mirror symmetries (t, x, y, ipi,ip2) (—t,—x,y,ij)i,ip2), 
(t,x,y^i^ 2 ) •->■ (t,x,-y,-ipi,-i()2) and (t, x, y, ipi, fo) ^ (t,x,y,ip 2 ,ipi). Under the change of 
the "barotropic/baroclinic variables" (ip + ,ip~) by the "layered variables" (tp 1 ,^ 2 ), these discrete 
symmetries are changed to the transformations mapping (i, x, y, to (— t, —x, y, 

(t, x, —y, —ip + , —ip~) and (i, x, y, —ip~), respectively. They exhaust the independent discrete 
symmetries of system (3) up to mutual composing and composing with continuous symmetries. 
Although the proof of this claim involves cumbersome calculations, the discrete symmetries are 
not essential for further consideration and hence will be neglected here. At the same time, they 
may be important for other purposes [3, 14]. 

5 Inequivalent subalgebras 

To classify the inequivalent subalgebras of the maximal Lie invariance algebra g of the system (1) 
(resp. the system (3)), we need the adjoint representation of the corresponding Lie group on g. 
See, e.g., [9, 21, 22] for details of classification techniques. We list only the nontrivial actions 
associated with basis elements of g: 



The classification of inequivalent one-dimensional subalgebras is straightforward. An optimal 
set of such subalgebras reads 



A\ = (d t + ad v + KF), Al = (d y + X(f) + bT), A\ = (X(f) + Z{g) + bF), (5) 



where a, b G R and / and g are arbitrary functions of t. 

The classification of the two-dimensional subalgebras yields the following list: 



A\ = (d t + kF, d y + X(u) + Z(fi) + pF), 

A\ = (d t + ud y + kF, X(e at ) + Z(uate at )),a ^ 

A\ = (d t + ud y + KF,Z(e at )), 

A\ = (d t + udy + kF , X(l) + Z(fi) + P F), 

A 2 _ 2 = (d t + udy + kF, Z{1) + P F), A 2 _ 3 = (d t + ud y ,F), 

A\ = (d y + X(f) + kF , ^(1) + Z(g) + P F), k P = 0, 

A 2 _ 4 = (d y + X(f),Z(g) + F), A 2 _ 5 = (d y + X(f) + kF, Z(g)), 

A\ = (X(f) + Zig 1 ) + kF, X(f 2 ) + Z(g 2 ) + pF), 



where a, b, k, p, u, p, a = const and /, / , f 2 , g, g 1 and g 2 are arbitrary functions of t. In 
the final subalgebra the tuples (f \g ,n) and (f 2 ,g 2 ,p) must be linearly independent for the 
subalgebra to really be two-dimensional. By the same reason, the parameter function g is not 



In fact, the above subalgebras are not single subalgebras but rather represent parameterized 
classes of subalgebras. This is why it would be beneficial to indicate the list of parameters in 
the notation of the corresponding algebras, but for the sake of brevity we omit this whenever 
possible. For the classes A 2 , A 2 _ A , A 2 , A^ 5 and A 2 _§ the adjoint actions and linear combinations 
of basis elements induce equivalence relations on the corresponding sets of parameters. For 
example, the subalgebras A 2 _ A (f, g) and A 2 _ 4 (f,g) are equivalent if and only if f(t) = f(t — e) 
and g(t) = g(t — e) for some e € R. 



M(e eZ{9) )d t = d t +eZ(g'), 
Ad(e eX ^)d t = d t + eX(f'), 
Ad(e sX W)d y = d y - eZ(f') 



Ad(e^)Z(g) = Z(g(t-e)), 
Ad(e^)X(f) = X(f(t-e)), 
Ad(e £d y)X(f) = X(f) + eZ(f). 



identically equal to zero in J?_ 
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6 Invariant reduction with one-dimensional subalgebras 

We present the complete list of submodels obtained under reduction using the list (5). For 
each submodel, we again determine its Lie symmetries, thereby seeking for hidden symmetries 
of the initial model. For a general discussion of the problem of hidden symmetries, see e.g. [1]. 
Throughout this section v , v 2 , v + and v~ will be assumed as functions of p and q. 

6.1 Subalgebra A\ 

A suitable ansatz for reduction of (1) under the first subalgebra of (5) is given by if) 1 = v 1 + bt, 
i\) 2 = v 2 — bt, where p = x, q = y — at. The corresponding reduced equations read 

awl - Fa H - v D + 2Fh - v lH + p + Fv l) + v lH + Fv l) = °> 

aw 2 + Fa{v\ - v 2 q ) - 2Fb - v 2 (w 2 + f3 + Fv\) + v 2 (w 2 p + Fv\) = 0, 

where w l = v l pp + v qq , i = 1,2. The maximal Lie invariance algebra of this system is 0i = 
(d p ,d q ,d v i,d v 2)- All operators from the algebra gi are induced by Lie symmetry operators of 
the original system (1) and hence there are no purely hidden symmetries. This is why we do 
not have to further reduce the above system by using the Lie method. The Lie reductions 
of the reduced system with respect to one-dimensional subalgebras of Qi are equivalent to Lie 
reductions of system (1) with respect to one of the listed two-dimensional subalgebras of q. The 
two-dimensional reductions of system (1) are exhaustively discussed in Section 7. 

6.2 Subalgebra A\ 

Reduction using A\. An ansatz associated with this subalgebra reads ip l = v 1 — \f'y 2 + by 
and if) 2 = v 2 — \ j'y 2 — by, where p = x — f(t)y and q = t. It reduces system (1) to the system 

{Hvl p ) q - f" - F{v\ - v 2 ) - bF{vl + v 2 p ) - bHvl pp + = 0, 
(Hv 2 pp ) q - f" + F{v\ - v 2 ) + bF{vl + v 2 p ) + bHv 2 ppp + 0t£ = 0, 

where it was convenient to introduce the new notation 

H = l + f 2 . 

To simplify the reduced system, we use the above mentioned barotropic/baroclinic variables, 
which is particularly obvious for this submodel, since it is a linear system of differential equations. 
By introducing w = v 1 + v 2 and v = v 1 — v 2 we are able to rewrite the resulting system via: 

{Hw P p) q - 2f" - bHvp P p + f3wp = 0, ^ 
(Hv P p) q — 2Fv q — 2bFw p — 2bHw ppp + (3v p = 0. 

Note that system (7) may be derived directly by means of reduction of (3) under the ansatz 
ip+ = w(p, q) — f'y 2 and = v(p, q) + 2by, where p and q are defined as above. 

The resulting system is now simplified in a way similar as presented in [4]. Namely, we 
integrate once the first equation with respect to p yielding 

(Hw p )q - 2f"p - bHvpp + I3w + h{q) = 0, 

where h is an arbitrary function of q = t. Then we apply the transformations of the unknown 
functions 

- jHf")> , 2f"p h A 2b f 



8 



and obtain the following system: 

{Hw p ) q +/3w- bHv pp = 0, ^ 
{Hv pp ) q - 2Fv q + j3v p - 2b(Hwp P p + Fw p ) = 0. 

First of all, we determine the Lie symmetries of this system. As system (8) in fact is a class 
of systems parameterized by the arbitrary function / = f(q) and the arbitrary constant b, 
it is necessary to solve a group classification problem [22, 29, 30]. That is, for the complete 
description of Lie symmetries it is necessary to seek for possible extensions of the Lie invariance 
algebras for special values of the parameters / and b, respectively. Recall that (3 and F are 
constant parameters which can be set to 1, so it is not required to also take into account the 
classification problem with respect to f3 and F. 

Group classification of the reduced systems. Conventionally, the first step in the procedure 
of group classification is the identification of the kernel G ker of the maximal Lie invariance groups 
of systems from class (8), i.e. the group which is admitted for any value of / and b. The Lie 
algebra g kcr corresponding to G ker can be obtained by solving the determining equations for 
Lie symmetries under the assumption of arbitrariness of / and b. The part of the determining 
equations not including / and b can be immediately integrated yielding 

£ p = ap + c, £ q = aq + d, rf = k\V + zip, q), rf" = k 2 w + z 2 (p, q), 

which are the coefficients of the most general infinitesimal generator of Lie symmetries £ p d p + 
t, q d q + rfdi, + rfd^, where a, c, d, k\, k 2 = const. The part of the determining equations 
explicitly including the parameters of (8) (the classifying part) in turn is: 

(aq + d)H' - 2aH = 0, (aq + d)H" - aH' = 0, 

aHH' + (aq + d)HH" - {aq + d)H 12 = 0, b(h - k 2 ) = 0, 

{Hz 2 p ) q + pz 2 - bHzlp = 0, 

{Hz l pp ) q - 2Fz] + /3z$ - 2b(Fz 2 + Hz 2 ppp ) = 0. 

It is straightforward to recover system (8) in the two last equations of system (9). For the 
general values of / and b splitting of the above system yields the conditions a = 0, d = and 
k\ = &2 and hence gives rise to the Lie invariance algebra generated by the operators 

dp, 1 = vdi + wdfi, C{z l ,z 2 ) = z 1 (p,q)d i) + z 2 (p,q)du,, 

where functions z 1 and z 2 run through the set of solutions of the system (8) for the fixed values / 
and b. The Lie symmetry operators I and Ciz 1 , z 2 ) arise due to linearity of (8). 

To investigate the problem of induced symmetries of system (8), we need to consider the 
same problem for system (7) at first. Up to linear combining, for general values of / and b the 
Lie symmetry operators of (7) induced by operators from g are exhausted by the operators d p , 
2dv and g(q)d$, which are induced by ^(1), F and 2(g), respectively. Here g runs through 
the set of smooth functions of q. It is obvious that the operator d y + X(f) + bT induces the 
zero operator. Additionally, if / = const the operator dt induces d q . Under integration of the 
first equation of system (7), Lie symmetry transformations generated by g(q)du, for any fixed g 
become equivalence transformations of the resulting system. By the above transformation to the 
unknown functions v and w, we gauge the function h arising under integration to zero and hence 
break the invariance of system (8) with respect to operators of the form g(q)dw- This is why 
any operator from 6 lying in the complement of the linear span of the operators d p , 2dy and 
g{q)dw and, additionally, d q in the case of / = const is a hidden symmetry of the initial system. 

The kernel algebra cj ker of class (8) is generated only by the operators d p , I and £(1, 0). This 
is because the set of generators Ciz 1 , z 2 ) is different for every representative of the class (8) as the 
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form of systems from the class depends on values of / and 6. However, it is not feasible to linearly 
combine solutions of different systems from the class (8) that contradicts belonging of /^(z 1 , z 2 ) 
to g ker for arbitrary values of z 1 and z 2 . Therefore, under solving the group classification problem 
for the class (8) it is natural to investigate extensions with respect to g| e ^ rather than with 
respect to g kcr . In other words, we should find all inequivalent values of the arbitrary elements / 
and 6 for which g^ is not the maximal Lie invariance algebra of the corresponding system of 
the form (8). Here the inequivalence is to be understood with respect to the equivalence group 
of the class (8). 

For this purpose, it is necessary to solve the classifying part (9) of the determining equations 
by taking into account for which forms of / and b an extension of g^ is admitted. It is obvious 
that for 6 = 0, the generator X from g^ 1 splits into the two generators vd$ and wd w . This 
splitting of X corresponds to the decoupling of the two equations (8). As the remaining classifying 
part of the determining equations is independent of 6, the extensions possible for different values 
of H are essentially not affected whether or not the two equations (8) are coupled. In case 
of b = we simply consider extensions with respect to q^ C q = (d p , vd$, wd w , ^-{z 1 , z 2 )) rather 
than to Q g fb' ^ * s cruc ial to remark that from the first three equations of system (9) only the 
first equation is independent. The other two equations are its differential consequences. As the 
investigation of extensions must be done up to equivalence, it would be necessary to compute 
the equivalence group of the class (8). However, it is obvious that this class admits scalings and 
shifts of q as equivalence transformations. For this reason, we only have to distinguish between 
the cases a / and a = 0. In the case a / 0, we can set a = 1 and d = by dividing the 
equation by a and shifting of q. As a result, we have H = xq 2 , where x is a positive constant 
in view of the definition of H, i.e., / = i-^/xc/ 2 — 1. The extension of g^ 1 is then given by the 
basis element pd p + qd q which is a hidden symmetry of the initial system (3). If a = 0, we have 
H = const and, moreover, H ^ 1, i.e., / = const and g™^ x = g^b + R- ecan that operator 
d q in this case is induced by the operator dt- Formally, we can scale the constant value H by a 
uniform scaling in p and q but this is not related to the nature of the problem. 

The uncoupled system (b = 0). We now proceed by studying the case 6 = 0, which leads to 
a decoupling of system (8): 

{Hw p ) q + /3u> = 0, (Hv pp ) q - 2Fv q + Pv p = 0. (10) 

The change of variables 

— -r, v = H(q)v, w = H(q)w, 

allows to transform this system to 

w pq + Pw = 0, v m -2F(H- 1 v)q + l3v p = 0. (11) 

Analogously to [4], in the first equation we recover the Klein-Gordon equation in light-cone 
variables. For this system, we only have one possibility for Lie reduction in the general case of H. 
Namely, we can reduce (10) by using the subalgebra (dp + \\vdy + \2wdw), where we take into 
account that for the decoupled case 6 = the generator vd% + wd w of (8) splits into the two 
single operators vdy and wd w (note that vdy = vdy and wd w = wd w ). The ansatz for reduction 
reads v = v(q)e Xl ^ and w = w(q)e X2 ^. It leads to the following system of ordinary differential 
equations 

X 2 w q + /3w = 0, \\v q - 2F(H~ l v) q + j3X x v = 0. 

The integration of the resulting system, the substitution of the obtained solution to the ansatz 
and the inverse change of variables yields 

. ci A P fdq\ A c 2 / [2H q F + p\ 1 Hdq 

w = — exp X 2 p — — / — , v = — exp Xip + ' 



H """" V w \2JHJ 1 H ~' F V J 2F- \ 2 H H 
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Recall that H = 1 + f 2 - For the values / = —l/k, where the constants k and I are wave 
numbers and Ai = A2 = ik, i 2 = — 1, this solution reduces to the well-known Rossby waves 
in the two-layer model. The solution for the barotropic mode w describes a single barotropic 
Rossby wave, which is independent of the vertical structure of the two-layer setting. In turn, 
the solution for v describes the evolution of the first baroclinic mode and explicitly depends on 
the vertical layer structure via the parameter dependency on F, which accounts for the density 
difference between the layers. Note that for a more general choice of the function /, this solution 
allows to derive series of wave solutions in a similar way as it was possible for the barotropic 
vorticity equation [4]. However, for / ^ const there arises an additional term proportional to y 2 
in the solution of ip + , which may violate the boundary conditions. Although such a violation of 
the boundaries does not exist for the solution of ip~ , the global (i.e. large-scale) realization of 
generalized Rossby waves is at once limited due to this restriction. 

The case b = shows that the classical Rossby wave solution can be recovered in two steps: 
Firstly reducing with respect to the operator d y + X(f) + bT and secondly performing reduction 
to a system of ordinary differential equations using a hidden symmetry of the submodel received 
in the first step. 

Reduction with respect to the additional symmetries for special values f(q) = const and 
f(q) = ±-y/ xq 2 — 1 will not be considered here, because for a decoupled system it is be better to 
construct exact solutions of each equation separately and then compose them to a solution of the 
entire system. In our particular case, the single Klein-Gordon equation has a wider maximal Lie 
invariance algebra than system (11), given by (dp, dq, pdp — qdq, wdw, z(p, q)dw), where g runs 
through the set of solutions of the Klein-Gordon equation. That is, we have more possibilities 
for finding exact solutions by Lie methods. Fortunately, it is not necessary to do this in view of 
large classes of exact solutions already known for the Klein-Gordon equation [25]. 

For the split system (11), it remains to determine the Lie symmetries and perform Lie reduc- 
tions of the second equation, 



where A = H 1 > 0. The determining equations for the coefficients of the Lie symmetry 
operator Q = £pdp + £pdp + rjdy of Eq. (12) not involving A can be integrated to give 

£ p = — ap + c, £ 9 = aq + d, rj = kv + z(p, q). 

The remaining classifying part of the determining equations reads 



Again, the second equation is a differential consequence of the first equation. For general A, 
splitting of (13) leads to the two essential Lie symmetry generators dp and vdy together with 
the linearity operators C(z) = z(p,q)dv, where z runs through the set of solutions of (12). 

Upon looking for Lie symmetry extensions, we should consider the cases where (a, d) 7^ (0, 0). 
We then distinguish two cases, (i) a ^ 0. We scale a = 1 and shift q to set d = and obtain the 
value A = xq 2 with the additional Lie symmetry generator qdq — pdp, where k is a nonvanishing 
constant, (ii) a = 0, d 7^ 0. In this case we find A = const and the additional Lie symmetry 
generator dq. The first operator again is a hidden symmetry, while the second generator is 
induced by dt- Note that the case of A = would yield wider symmetry extensions but cannot be 
realized in the present case since by definition A ^ 0. Both the cases of extensions correspond to 
those above for arbitrary b by taking into account the change of variables leading to system (11). 

We now consider the Lie reductions of Eq. (12). For general A, the only nontrivial possibility 
for reduction is given by dp + Xvd^, which was already considered above. The solution is 





(aq + d)Aq - 2aA = 0, aAAq - (aq + d)A 2 ^ + (aq + d)AA q - q - = 0. 



(13) 
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which can be combined with arbitrary solutions of the Klein-Gordon equation to yield a solution 
of the decoupled system (11). It remains to study Lie reductions involving symmetry extensions. 

In the first case we have A = xq 2 . The maximal Lie invariance algebra of the equation (12) 
with A = xq 2 is given by (qd q — pdp, dp, vd v , C{z)). One more inequivalent one-dimensional 
subalgebra of this algebra is appropriate for Lie reduction. It reads (qd q — pdp + Xvdy). An 
ansatz for the reduction is v = v(r)q x , where r = pq. The corresponding reduced equation is 

rv rrr + (A + 2)v rr + (/3 - 2xFr)v r - 2xF(X + 2)v = 0. 

For A = —2, this equation is solved in terms of Whittaker functions M mjTl (r), W m ^ n (r): 

v = dj M m , n (V8^Fr) dr + C 2 J W m , n {V^Fr) dr + C 3 , 

where m = (5{d>xF)~ 1 / 2 and n = 1/2. Moreover, for A = —k, k £ N, this equation admits 
polynomial solutions. 

In the second case of extension, we have A = const. Then, the maximal Lie invariance algebra 
coincides with {d q ,dp,vdy,C(z)). Again, one more inequivalent one-dimensional subalgebra can 
be used to carry out Lie reduction, which is (d q + ndp + Xvdv). An appropriate ansatz for 
reduction is v = v(r)e Xq , where r = p — nq. Plugging this ansatz into (12), we find 

KVrrr ~ Xv rr ~ (2AFk + (3)v r + 2AF\t) = 0. 

This is a linear third-order ordinary differential equation with constant coefficients and thus 
can be solved by standard methods. In particular, as k, A and A are arbitrary constants, we 
can determine them upon prescribing a solution of the associated characteristic equation. This 
allows one to generate wide classes of solutions with rather different type, such as e.g. periodic 
wave solutions. 

Since (12) is a linear partial differential equation, using the same symmetries and the method 
of extended Lie reduction described in Appendix A we can obtain much wider families of exact 
solutions for this equation, cf. the example in this appendix. 

The coupled system (b ^ 0). For the coupled case b ^ the Lie reduction of system (8) is 
quite similar to the case 6 = 0. For the sake of completeness, we list here all the reduced models 
that are possible for different values of H. 

For general H, the only possibility for reduction is due to the generator Q = d p + XT. 
The corresponding reduction ansatz is v = v(q)e Xp , w = w{q)e' Kp . Plugging the ansatz into 
system (8), one obtains 

XHw q + (XH q + p)w - bX 2 Hv = 0, 

{X 2 H - 2F)v q + (X 2 H q + f3X)v - 2bX{X 2 + F)w = 0. 

For H = xq 2 , the additional inequivalent reduction using Q = pd p + qd q + XX is possible. 
Utilizing the ansatz v = v(r)q x , w = w(r)q x , where r = pq~ l , leads to the system 

k(X + l)w r — xrw rr + j3w — bxv rr = 0, 

xrv rrr — X>cv rr + 2F(Xv — rv r ) — j3v r + 2b(xw rrr + Fw r ) = 0. 

For H = const we can also reduce (8) using Q = d q + k8 p + XX. The ansatz for reduction is 
v = v(r)e Xq , w = w{r)e Xq , where r = p — nq. The resulting submodel is 

H{nw Tr — Xw r ) — (3w + bHv rr = 0, 

H(nv r rr — Xv rr ) — (2Fk + 0)v r + 2FXv + 2b(Hw r rr + Fw r ) = 0. 

All of the above reduced systems are linear systems of ordinary differential equations and, 
moreover, the last system has constant coefficients. As in the decoupled case 6 = 0, for finding 
more exact solutions we could also apply the method of extended Lie reduction to (8) with 6^0. 
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6.3 Subalgebra A\ 

Reduction using ^,3. For this subalgebra, it is convenient to start with barotropic/baroclinic 
variables from the beginning. Since in the case / = no Lie reduction is possible, we assume 
that / 7^ 0. An appropriate ansatz for reduction then reads 

ip + = v + — 2 ^ ^ — -x, ip~ = v~ + ^yX, 

where p = y and q = t. Plugging this ansatz into (3) gives reduction to the system: 

+ / P ~ 9 , + , 2 8) + -v- =0 
ppg j \ u ppp 1 ^HJ 1 j ppp 1 



V 



II 1 

2Fv- - J -^j^- (v m - 2Fv~ ) + - (v+p + 2Fv^ + 2/3) = 0, (14a) 



where it can be seen that the coupling between the barotropic and baroclinic parts is again 
provided only due to the existence of generator T . To solve this system, we integrate the first 
equation twice with respect to p to yield 

v+ ~ ! -^J JL 4 + 2jv+ + h - f v- - \ f j^ + P -fp 2 + h\q)p + h°(q) = 0, (14b) 

where h 1 and h° are arbitrary smooth functions of q. By means of the change of unknown 
functions 



v 

where 

2 



y+ + 7 2 p 3 + 7 i ^ p + 7 ° ( g ) ; v~ =v~ +5 2 (q)p 2 + 6 1 (q)p + 8°(q), 



l l = ~jj (2^ 2 + fh l )dq, 5 1 = -2f J ^dq, 



1 r 1 1 n 1 r hF i l - 9 F ^ + ffi 

7 7 s J f{91 + + f )dq ' fJ — 1 

we are able to reduce system (14) to the corresponding homogeneous form: 

Vq J Vp + 2jV + + jV p = 0, 



Vp Pq ~ 2Fv~ — (Vp pp - 2Fv~) + j(Vp PP + 2Fv+) = 0. 

System (15) is simplified further using the transformation 

p = f(q)p- J g(q)dq, q = q, v + = v + , v~ =v 

In the new variables, system (15) becomes 

(f 2 V + )q+b(f 2 V~)p = 0, 

(f%~ - 2Fv~)q + b(f 2 v± + 2Fv+)p = 0. 

The first equation of (16) can be used to introduce a potential variable, via Vp = f 2 v + and 
Vq = —bf 2 v~. Upon introducing q = f f 2 dq, the second equation then becomes 

c2/j-2ta.. \ oz7^2 T / f.2 ' — 2F ^ 



(15) 



(16) 



f 2 (fVppq)q ~ 2Ff%q - b 2 ( Vpppp + j^Vpp ) = 0. (17) 
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The decoupled system (b = 0). The general solution of the decoupled system (16) is 



r 

where C 1 and C 2 are arbitrary functions of p = f(q)p — f g{q)dq and i? 1 and # 2 are arbitrary 
functions of q. 

The coupled system (6^0). As we found the general solution for the decoupled case of 
systems of the form (16), we should not study further Lie reductions for this case. Hence we 
carry out Lie reductions solely for the case of b ^ 0. For this reason, it is necessary to solve 
the group classification problem for the class (16) with /, g and b as arbitrary elements under 
the assumption that 6/0. As the class (16) consists of linear systems, the procedure of its 
group classification should be done as described in the previous section. The general form of 
Lie symmetry operators of a system from the class (16) is £, p dp + i q d q + r] + d^+ + T]~d%- , where 
the coefficients are functions of p, q, v + and v~ . The solution of the corresponding determining 
equations gives 

£ p = ap + c, £ 9 = aq + d, T] + = kv + + z(q, p), rf~ = kv~ + z 2 (q,p). 

Additionally there is the single classifying equation (aq + d)f q — af = (up to its differential 
consequences). For the general value of the arbitrary element / we obtain a = d = 0. The 
maximal Lie invariance algebra in this case reads 

Of" = (dp, l = v + d i)+ +d~dv-, C(z l ,z 2 ) = z 1 (p,q)dj,+ + z 2 (p, q)d$- ), 

where (z 1 ,^ 2 ) run through the set of solutions of system (16). Inequivalent extensions of the 
algebra jjy en only exist for / = const and / = x.q. The first case leads to the extension of 0j Cn 
by dg, the second case gives the extension of gj Cn by pdp + qdq. 

To study the problem of induced Lie symmetries of system (16) with b / 0, we need to 
consider the same problem for the initial reduced system at first. Up to linear combining, for 
general values of / and g the Lie symmetry operators of (7) induced by operators from g are 
exhausted by the operators 

fq 2„ , 9 A Q n b 



*(/) = -2 if q p ~ Jfp + jf) d v+ - 2-fd v - , F=2d v -, Z(~g) = ~g(q)d v+ , 

which are induced by <-f (/), T and Z(g), respectively. Here / and g run through the set of smooth 
functions of q = t. Additionally, if / = const the operator d y induces d p and if also g = const 
the operator dt induces d q . Lie symmetry transformations generated by X(f) + Z (g) with any 
fixed / and g are equivalence transformations of the system (14). By these transformations, we 
gauge the functions h° and h 1 in (14b) to zero and hence break the invariance with respect to 
the operators X (f) and Z(g) with the general values of / and g. Summing up, we can say that 
almost all operators from g^ 011 and its extensions are hidden symmetries of the initial system (1). 

It was shown before that by introducing a potential variable, a system of the form (16) can 
be converted into a single equation in the conserved form (17). Since such an equation may have 
additional symmetries compared to the original system, we also carry out the group classification 
of the class (17) (again under the assumption that b ^ 0). Solving the determining equations 
for the coefficients of a Lie symmetry operator ^dp + ^dq + r) V dy of (17) gives 

£P = ap + c, ^ = 3aq + d, 7] V = aV + z(p, q), 

together with the classifying equation (3aq + d)f — af = 0. For arbitrary /, this imme- 
diately implies that a = d = and hence gives rise to the maximal Lie invariance algebra 
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{dp, Vdy, z{p,q)dy), where g is an arbitrary solution of (17). There are two inequivalent ex- 
tensions of this algebra, depending on either a ^ 0, d = or a = 0, d ^ 0. Since scalings and 
shifts in q are equivalence transformations, we may first set a = 1/3, d = with / = x^fq. The 
additional generator then reads pdp + 3qdq. The second case of extension is given by a = and 
d = 1, leading to / = const with the corresponding additional generator dq. 

Comparing this result with the classification of system (16), we see that all Lie symmetries 
of the potential equation (17) are induced by Lie symmetries of (16) (note the change of the 
variable q in the potential case). That is, for the reduced system (16) no purely potential 
symmetries associated with the potential equation (17) exist. 

Now that we have investigated all symmetry extensions for particular values of f{q), it remains 
to present the corresponding Lie reductions of (16). The only feasible way for reduction for 
general / is due to the operator dp + XI. The ansatz for reduction then is v + = v + [q)e x P and 
v~ = v~(q)e x P. The system of reduced equation reads 



bf 2 Xv~ = 0, (f 2 X 2 v~ - 2Fv~ )g + bX ( A 2 + ^ ) 



where u = f 2 v + . The case of A = gives only a trivial solution and will not be considered here. 
For A it is possible to solve the first equation for v~ . Plugging the corresponding expression 
into the second equation, the following homogeneous second order ordinary differential equation 
with variable coefficients for u is obtained 



4i? ll u ^ 2 A 2 / 2 + 2F 
X 2 P-2F f 9 X 2 f 2 -2F 



Two more Lie reductions are possible for the particular values / = n = const and / = xq. 
In the first case, the maximal Lie invariance algebra reads {dq, dp, I, C^ 1 , z 2 )). We aim to 
reduce (16) using the generator Q = dq + ndp + XT. The ansatz for reduction is v + = v + {r)e x ^ 
and v~ = v~(r)e x ^, where r = p — nq. Plugging this ansatz into (16), we obtain 

kv^~ — Xv + — bv~ = 0, 

x 2 {nv~ rr — Xv~ r ) — 2F(kv~ — Xv~) — b{x 2 v^ rr + 2Fv^) = 0. 

This is a coupled system of first and third order ordinary differential equations with constant 
coefficients and thus can be solved explicitly using standard techniques. 

In the second case, the maximal Lie invariance algebra is given by {qdq+pdp, dp, I, C{z l , z 2 )). 
One-dimensional reduction is feasible using the generator Q = qdq +pdp + XI. The corresponding 
ansatz for reduction is v + = v + (r)q x and v~ = v~(r)q x , where r = pq" 1 . This ansatz leads to 
the following reduction of system (16): 

ru+ - (A + 2)v + - bv~ = 0, 

rv~. r - (2x 2 + x(A - 2))u~ - 2F{rv~ - Xv~) - b(x 2 v+ r + 2Fv+) = 0. 

Remark 3. Substituting back the solutions of the reduced equations into the ansatz for the 
original unknown functions ifj + and ifj~ , we find that all solutions have time-dependent polyno- 
mial parts in x and y. Since this part is not compatible with typical boundaries of the two-layer 
equations (as discussed in Section 8), only solutions for the restricted case of b = g = 0, / = const 
might give candidate solutions that could be realized in the framework of geophysical fluid dy- 
namics. However, since the solution of the decoupled case 6 = for g = 0, / = const is rather 
trivial, it is not too interesting from the physical point of view. 
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7 Invariant reduction with two-dimensional subalgebras 



Now we carry out Lie reductions of (1) or (3) to systems of ordinary differential equations using 
the optimal set of two-dimensional inequivalent subalgebras. Note that not all subalgebras give 
rise to classical group-invariant reductions. Namely, all the subalgebras from the list (6) with 
negative subscripts cannot be used for the construction of well-defined ansatzes for the dependent 
variables. However, these algebras would be well suited for the construction of partially-invariant 
solutions, but we do not pursue this idea further in the present paper. 

7.1 Subalgebra A\ 

An ansatz for group-invariant reduction using this subalgebra is ip l = v 1 (p) + nt + {p + p)y and 
■tp 2 = v 2 {p) — Kt + (p — p)y, where p = x — vy. System (1) is reduced by this ansatz to 

-(p + fi)(l + v 2 )v\ pp + Fp(Vp - v 2 p ) - Fp(vl + v 2 p ) - 2Fk + ^ = 0, 
(p - + v 2 )v 2 ppp - Fp(y x p - v 2 p ) + Fpivl + v 2 p ) + 2Fk + f3v 2 = 0. 

Integrating once with respect to p, the above system becomes an inhomogeneous system of two 
second order linear ordinary differential equations with constants coefficients, provided we screen 
out the singular cases oi p = p = 0, p = p > 0, p = p < 0, p = —p > and p = —p < 0. 
The solution of this system in the nonsingular case is straightforward but a bit lengthy and is 
therefore omitted here. Hence, we will focus on the listed singular cases only. 

Case p = p = 0. The general solution for this case is 

t 2Fk 2 2Fk 
v = — j-p + ci, v = — p + c 2 , 

P P 

which in the original variables (tp , tp ) gives a solution linear in (x,y) and thus represents a 
constant wind field in both layers. 

Case p = p > 0. This case leads to the semi-coupled system 

2p(l + v 2 )v\ pp + 2Fpv 2 p + 2Fk - pvl = 0, 2Fpv 2 p + 2Fk + f3v 2 p = 0, 
which is integrated to yield 



yl = C1 exp (v mtW) + C2 exp (" V Mr^-f) + + C3 ' 

2 2Fk P 

In the original variables, this represents a simple exponential solution and is thus unphysical. 
Case p = p < 0. The general solution in this case is 



v 



2 2Fnp 
V = 2-FjpJ^ + C ^ 

which in the original variables represents a single (stationary) Rossby-wave in the upper layer 
and a constant velocity field in the lower layer, see Figure 1. This is a typical situation, in the 
study of baroclinic instability: An initial disturbance in the middle of the troposphere may start 
to grow while the lower part of the troposphere does not exhibit any peculiarities. It is not 
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before the upper Rossby-wave starts unstable growth, that the lower layer also begins to show 
some wave-like disturbances, which subsequently may lead to the onset of cyclogenesis. Hence, 
the above exact solution may characterize the situation at the onset of baroclinic instability, 
where due to external forcing a Rossby wave in the upper layer is generated, while the wind 
field in the lower layer is still unaffected. 




Figure 1: Rossby wave solution in the quasi-geostrophic two-layer model. 

The cases p = —fx > and p = —fi < give the same solutions as in the two previous cases, 
except for interchanging the two layers, i.e. v 1 <->■ v 2 and permuting the sign of the linear in p 
term. 

7.2 Subalgebra A\ 

It is convenient to use at once the quasi-geostrophic equations in terms of barotropic/baroclinic 
variables to perform the reduction. The ansatz we choose is: ip + = v 1 (p) — 2apx and = 
v 2 (j>) + 2nt, where p = y — vt. Using these variables, the resulting submodel of (1) is: 

{v + ap)v l ppp + 2b/3p = 0, (v + ap)v 2 ppp - 2F{y + ap)v 2 + 4Fk = 0, 

which is a decoupled system of third order linear ODEs. The case of b = is trivial and will 
not be considered here. For 5^0 the general solution is 

i v n , , . / 9 v v 2 \ 1 „ / o 9is 9 6v 2 3v 3 \ 

V = HV + ap) ( P +2 a P+ ^)~f [ P + Ya P + ^ 2 - p+ 2^) 
+ cip 2 + c 2 p + c 3 , 

v 2 = ^ (2hi(u + ap) + Ei + ap)) e^ v+ ^ + Ei (-^(i/ + ap)) e"^^)) 

+ c 4 e^ + C5 e-^ + C6 , 

where Ei(z) = t~ 1 e~ t dt denotes the exponential integral. In terms of the ip + /?p~ variables 
this solution is the superposition of some polynomial with an exponential function in y-direction. 
From the meteorological point of view, this solution does not seem to be relevant. 

7.3 Subalgebra A\ 

Using the barotropic/baroclinic variables, the ansatz for reduction under this subalgebra is 
ip + = v 1 (p) + 2px, ip~ = v 2 {p) + 2nt + 2px, where p = y — vt. The reduced system then is 

- V) v \ vv ~ P v l PP ~ 2 PP = 0, 

(y - p)v 2 pp - pv l ppp - 2F(v - p)v 2 p - 2Fpv\ + AFk - 2fip = 0. 
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This system is now completely integrable. The case p = and /i / 1/ leads to a decoupled 
system of equations which can be integrated easily. The same also holds in the case p 7^ and 
v = p. Hence, we focus on the case where p ^ and v and p are arbitrary. The first equation 
can be integrated at once three times, yielding the relation between v 1 and v 2 , given by 



v 2 



i f(p - pjv 1 - ^(3pp 3 + cip 2 + c 2 p + c 3 ^j 



where ci, C2 and C3 are arbitrary constants. Then, integrating once the second equation and 
substituting the expression for v 2 produces a second order ordinary differential equation with 
constant coefficients, from which we determine v : 



v 1 



A 1 sm \l -Ijp + A 2 cos \l '-^p- ^ (5 3 p 3 + 5 2 p 2 + 5 1 p + 5°) + -^2 (65 3 p + 25 2 ) , 



where 



^ = _2F (-{v - p) 2 + pY 7 i = i( I/ _ Ai )2_ p 

\P ) P 

,3 2F PP, x r2 2 ^F. 

^ = -2f ^ + C2i? (,- M )-(2F.-/3 P )y jo = 2(cx-C3F) ^ 

VP / P 

provided that 7 2 /7 1 > 0. In the particular case of v = p, which leads to a considerable 
simplification of the above solution, this condition is verified. In the case of 7V7 1 < we 
can find a solution in terms of exponential functions, which is not presented here. Plugging 
this solution into the ansatz for the original unknown functions, the above solution gives the 
combination of a traveling wave in y-direction with a third order time-dependent polynomial in 
x and y. For usual fixed boundaries in north-south direction, this is an unphysical solution. 

7.4 Subalgebra A\ 

An appropriate ansatz for this subalgebra is 

ip + = v 1 {p) - f'y 2 -2g(fy - x), = v 2 (p) + 2ny - 2p(fy - x), p = t, 

where we have again employed the barotropic/baroclinic variables. This ansatz enables reduction 
of (3) to the system 

f-fa = 0, v 2 + 2 K g-^ = 0. 
The general solution of this system is 

v =6(p), v =-—f+—p + c, 

where 9 is an arbitrary function of p and c is an arbitrary constant. The first equation of 
the reduced system is the compatibility condition of the initial system (3) and the invariant 
surface condition corresponding to the subalgebra A 2 . It implies the constraint g = f"//3 for 
the parameter functions / and g, which is the necessary and sufficient condition for system (3) 
to have solutions invariant with respect to the subalgebra A 2 . This solution has no obvious 
physical importance in dynamic meteorology. The reason is that this is a simple polynomial 
solution with time-dependent coefficients. We note that the function 9 can be set to zero due to 
gauging of the stream functions generated by 2(g). 
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8 Invariant reduction of boundary value problems 



In this section, we aim to discuss admitted symmetries in the presence of boundaries. It is 
commonly assumed that group-invariant solutions may describe the behavior of a system that 
is far away from boundaries and hence a consideration of restrictions imposed by boundaries 
is usually omitted. However, as shown e.g. in [12, 18], there may be situations where the 
system without boundaries is not simply the limit of a system with very distant boundaries. 
Consequently, consideration of boundaries may be necessary even for a conceptual understanding 
of the model evolution. Moreover, as noted in the two previous sections, some of the group- 
invariant solutions corresponding to the optimal sets of inequivalent subalgebras give rise to 
unphysical solutions due to a violation of boundaries. We now compute those symmetries that 
are admitted by the boundaries and hence discuss which solutions could be compatible with the 
boundary value problem. 

In the atmospheric sciences, for equations on the /3-plane commonly a channel flow is assumed, 
which implies rigid boundaries y = and y = Y in north-south direction. In east-west direction, 
one usually assumes periodic boundaries at x = —L and x = L or an infinitely extended domain. 
In this setting, imposed conditions for the two-layer model are 



The second condition implies conservation of circulation at the boundaries. According to [7], 
for a boundary value problem to be invariant, three conditions must be satisfied: (i) invariance 
of the equation, (ii) invariance of the domain, (iii) invariance of the values on the boundaries. 
The first condition was already established in Section 3, so it remains to verify (ii) and (iii). 
Although it is possible to solve this problem on the stage of the Lie algebra using the infinitesimal 
method [7], we find it more comfortable to work with the finite group transformations. The most 
general form of a continuous symmetry transformation of (1) is given by 

(t,x,y,ipi,ip 2 ) >->■ (t + ei,x + f,y + e 2 ,ipi - f'y + g + e 3 , ^2 - f'y + g - £3), 

where ex, £2 and £3 are arbitrary constants and / and g are arbitrary smooth functions of t. 
This transformation has to preserve both the domain and the boundary values. We now discuss 
the channel flow with three possibilities for boundaries in east-west direction. 

Infinite domain. If there are no sidewalls in east-west direction, that is L — > 00, the most 
general symmetry preserving the boundary value problem is given by 

(t,x,y,tpi,tp2) ^ {t + ex,x + h,y,ipi - h'y + g + £3,^2 - h'y + g - e 3 ). 

where h is an arbitrary function of t with h" = 0. Hence, we have h = e^t + £5, where £3 
and £5 are arbitrary constants, leading to the important result that Galilean boosts preserve the 
boundary value problem. 

Periodic boundaries. For periodic boundaries, we have ipi(t, —L,y) = tpi(t,L,y). Similar 
calculations as above imply that the boundary-preserving symmetry group is the same as for an 
infinite-domain. 

This shows that the Rossby wave solution is admitted by the boundary value problem. This 
may serve as a "symmetry explanation" for the prominent occurrence of this solution in geo- 
physical fluid dynamics. 

Limited domain. The model of the limited domain in east-west direction is very natural 
in oceanography but can be also realized in the atmospheric sciences as flow in a mountainous 



L 




for 




(19) 
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region. For the purpose of simplicity, we assume a rectangular domain. Besides (19), this setting 
requires the additional conditions 

Y 



where L and Y denote the length and the width of the rectangle, respectively. Symmetries that 
are compatible with this boundary value problem are 

{t,x,y,ipi,ip 2 ) h> (t + e 1 ,x,y,ijji +g + e 3 ,^ 2 + g- e 3 ). 

Hence, the only group-invariant solution that can be realized on this domain is a stationary 
solution. 

9 Conclusion 

In this paper, we consider the baroclinic two-layer model from the viewpoint of Lie symmetries. 
The maximal Lie invariance algebra and the complete point symmetry group including both 
continuous and discrete symmetries are computed. We construct exhaustive sets of one- and 
two-dimensional subalgebras which are inequivalent with respect to the adjoint action and carry 
out the corresponding Lie reductions in one and two variables. This completely solves the 
classical problem of Lie reduction for the two-layer model. The procedure leads to various 
classes of exact solutions, some of which are well known in the atmospheric sciences, including 
barotropic and baroclinic Rossby waves. Finally, also the two-layer boundary value problem is 
investigated in the light of admitted Lie symmetries. We obtain a result, which is similar to that 
obtained in [32] (see also [22, pp. 379]) for the Navier-Stokes equations, namely that periodic 
boundary conditions admit Galilean boosts as symmetry transformations. 

Although there are still a large number of obviously unphysical solutions, the study of the 
classical Lie problem is a necessary first step for the consideration of partially invariant solutions 
and nonclassical symmetries, which we save for future investigations. In addition, these solutions 
are of undeniable value when it comes to a numerical implementation of the two-layer equations, 
which can employ several kinds of boundary conditions. The solutions constructed can be used as 
benchmark tests to assess the quality of the numerical scheme involved by addressing issues such 
as convergence rates and the reproduction of correct phase space velocities of wave-like solutions. 
Moreover, computing differential invariants of subalgebras of the maximal Lie invariance algebra 
determined in this paper, one can extend the set of exact solutions, e.g. by the construction of 
differentially invariant solutions [10]. 

From the mathematical point of view it is interesting to note the following properties of 
each of the above inequivalent Lie reductions, excluding that obtained using the subalgebra A\: 
The associated reduced system is linear and has a simpler form in the "barotropic/baroclinic 
variables". The coupling of the reduced equations written in terms of "barotropic/baroclinic 
variables" is only under the presence of the "baroclinic" operator T = d^i —d^2 in the subalgebra. 
The reduced equation for the baroclinic part of the model is structurally more complicated than 
that for the barotropic part. 

The linearity of the reduced system is especially beneficial in the second and third cases of 
reduction in one variable. The linear superposition principle available for the resulting systems of 
two-dimensional partial differential equations allows one to generate wide sets of exact solutions 
by linear combining of different solutions and substituting them back into the ansatz for the 
original unknown functions. Furthermore, there exist a number of special techniques for finding 
exact solutions of linear partial differential equations. One of them is presented in Appendix A. 

Since the two-layer model is only capable of resolving the barotropic mode and the first 
baroclinic mode, it would be interesting to study symmetry properties of multi-layer models. 
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This would be a preliminary step on the way to the investigation of a continuously stratified 
atmosphere. On the other hand, from the standpoint of application, a deeper investigation 
of layer models may even be more important than the three-dimensional system of governing 
equations. This is true since numerical utilization of these equations calls for some discretization, 
hence naturally leading back to the model of a layered atmosphere. Therefore, the present 
investigation of the two-layer model may not only be interesting for historical reasons. 
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A Extended Lie reduction of linear PDEs 

For clarity of the presentation we confine our consideration to the case of a single differential 
equation with a single dependent variable. 

Consider a linear partial differential equation C: Lu = in the unknown function u of 
n independent variables x = (xi, . . . ,x n ), where L is the associated linear differential operator. 
In what follows we use the summation convention for repeated indices. The indices i, j and k 
run from 1 to n, the indices a and b run from 1 to m. 

Suppose that the equation C possesses a nontrivial Lie symmetry operator Qq of the form 
Qo = C( x )d Xl + fl{x)ud u , where 7^ 0. Then for an arbitrary constant A the equation C 
obviously possesses also the vector field Q\ = Qq + Xud u as a nontrivial Lie symmetry operator. 
By Q\ we denote the differential operator which acts on functions of x and is associated with the 
operator Q\, i.e., Q\ = —^ t {x)d Xi + r](x) + A. For any m G N the differential function (Q\) m u 
is well known to be a characteristic of a generalized symmetry of C, and hence any associated 
generalized ansatz reduces the equation C to a system of m linear differential equations in m new 
unknown functions of n— 1 new independent variables invariant with respect to the operator Q\. 
In order to construct an ansatz, we should integrate the partial differential equation (Q\) m u = 0. 
The general solution of this equation gives the ansatz 

m am— a 

u = h(x)e xe y^ a (u J ) 7 ^ -, 

L^-r ' <m- a )\ 
o=l v ' 

where uj = (a; 1 (x) , . . . ,cj n ~ 1 (x)) is a tuple of functionally independent solutions of the equation 
C u Xi = 0, which are assumed to be invariant independent variables, 6 = 9{x) is a particular 
solution of the equation C u Xi = 1, h = h(x) is a particular nonvanishing solution of the equation 
C u xi = iju and <p a = </? a (w) play the role of new unknown functions. 

In view of the Lie invariance with respect to the operator Qq, the equation C is mapped by 
the point transformation 

xx=cj 1 (x), x n _i = uj n ' 1 (x), x n = 9(x), u = -^— 

n{x) 

to the equation Lu = 0, where the coefficients of the operator L do not depend on x n . In the 
new variables (x, u) the ansatz (20) takes the form 

u = e ^Y j v a {x l) ...,x n - l )-^— 

(m — an 

a=l v ' 

After substituting the ansatz (21) into the transformed equation Lu = 0, dividing the resulting 
equation by e Xxn and subsequently splitting with respect to different powers of the variable ot n , 
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(20) 



(21) 



we obtain, at least for the general value of A, the system 1Z of m differential equations with 
respect to the functions <p a in n — 1 independent variables (x\, . . . ,x n -i). In singular cases, for 
certain values of A some of the equations are identities. The same reduction is obtained by the 
substitution of the ansatz (20) into the initial equation C 

If the basic field is real, we can consider complex values of A, construct the corresponding 
complex exact solution and then take its real and imagine parts in order to obtain real solutions. 

The above consideration has a nice interpretation within the framework of Lie symmetries. 
Introducing the new dependent variables v a = (Q\) m ~ a u, instead of the single mth order linear 
partial differential equation (Q\) rn u = for finding a generalized ansatz, we obtain the system 
of m first order linear partial differential equations 

Q x v 1 = 0, Q x v a = v a -\ a = 2,...,m. 

As Q\ is a Lie symmetry operator of the equation C, each function v a satisfies this equation. 
To give the interpretation, we consider the system S of m copies of the initial equation C 

Lv 1 = 0, . . . , Lv m = 0. 

This system obviously possesses the operators Qq = ^ l {x)d Xi + rj(x)v a d v a and v b d v a as its Lie 
symmetry operators. Consider a linear combination of these operators, Qa = Qo + A.abV b d v °., 
which is also a Lie symmetry operator of the system S. Here and in what follows A a ;, are 
constants. Up to the equivalence generated by adjoint action of the Lie symmetry group of S on 
the corresponding Lie invariance algebra and due to the linear superposition principle, we can 
assume without loss of generality that the matrix A = (A a &) is the single m x m Jordan block 
with an eigenvalue A, 



A = J? = 
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The invariant surface condition for the operator Q\ with A = J7 1 consists of the equations 

?v 1 = (ri + \)v 1 , Cv a = (t] + X)v a + v a -\ a = 2,...,m, 
and an ansatz constructed with this operator has the form 

/ = %)e M ^V)7^, (22) 
6=1 ^ '' 

where the notation from the ansatz (21) is used. According to the general theory of Lie re- 
duction [21], the ansatz (22) necessarily reduces the system S to a system in the functions tp a , 
which obviously coincides with the system 1Z obtained by reducing the single equation C using 
the generalized ansatz (20). 

If the equation C is considered over the real field and the eigenvalue A is complex, it is not 
necessary to pass to the associated real Jordan block. As above, we can find the corresponding 
complex exact solution and then take its real and imagine parts in order to construct real 
solutions. This additionally justifies the usage of complex values of A in the real case. 

Example. For the general value of A, the second equation of system (11) admits only one inde- 
pendent nontrivial Lie symmetry operator, d p . Consider a system of m copies of this equation: 

v; pq -2(Av a ) g + v a p =0, (23) 
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where for simplicity we have omitted bars over the variables and scaled F = 1 and j3 = 1. This 
system admits the Lie symmetry operator Qa = d p + A a j,v b d v a ) where A = JP. The invariant 
surface condition associated with Qa reads 

uJ = At; 1 > i$ = \v a + v*~ 1 , a = 2, ...,m. 

Its general solution provides us with an appropriate ansatz for Lie reduction: 

„a— b 



7 a = exp(Ap) £y(. 



ip k = Cke'^ + e'^ I — fJ — dq, k = 3, ...,m, 



(a-b)V 
6=1 v ; 

where w = q = t is the invariant independent variable. Substituting this ansatz into system (23) 
yields the system of ordinary differential equations for ip a 

L^p 1 = 0, 

Ly? 2 + 2X^1 + ip 1 = 0, 

L/ + 2A(^~ 1 + p*" 1 + ^ 2 = 0, fc = 3, ... ,m, 

where the operator L is given by L := (A 2 — 2^4)9 g — 2A q + A. 
The solution of the above system is: 

ip 1 = cie _c , 

2 -c , -c /" ^ + 2A ^ C j 
V? = c 2 e <■ + e W 2A _ X 2 e d( l' 

2) 

2A-X 2 
where 

y 2A g - a 

C = J 2A^ dq - 

In the special cases A = const and A = x(q + Co) 2 , where x, Co = const, we can make more 
generalized reductions of the equation under consideration, which involves operators extending 
the Lie invariance algebra of the general case, cf. Section 6.2. 
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